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Abstract 

In this article, we study superfluid behavior of a gas of spatially indirect magnetoexcitons with reference to a 
system of two graphene layers embedded in a multilayer dielectric structure. The system is considered as an 
alternative of a double quantum well in a GaAs heterostructure. We determine a range of parameters (interlayer 
distance, dielectric constant, magnetic field, and gate voltage) where magnetoexciton superfluidity can be 
achieved. Temperature of superfluid transition is computed. A reduction of critical parameters caused by impurities 
is evaluated and critical impurity concentration is determined. 
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1 Introduction 

Recent progress in creation of heterostructures with two 
graphene layers separated by a thin dielectrics [1] opens 
possibilities to use graphene for creation of multiple 
quantum well structures with separately accessed con- 
ducting layers. In [1], Si0 2 substrate and A1 2 0 3 internal 
dielectric layer were used. Another promising dielectric 
is hexagonal BN [2], It has a number of advantages, 
such as an atomically smooth surface that is free of 
dangling bonds and charge traps, a lattice constant 
similar to that of graphite, and a large electronic 
bandgap. 

The attention to graphene heterostructures is caused, 
in some part, by the idea to use them for a realization 
of superfluidity of spatially indirect excitons [3-9]. 
Bound electron-hole pairs cannot carry electrical charge, 
but in bilayers they can provide a flow of oppositely 
directed electrical currents. Therefore, exciton super- 
fluidity in bilayers should manifest itself as a special 
kind of superconductivity-the counterflow one, that 
means infinite conductance under a flow of equal in 
modulus and oppositely directed currents in the layers. 

The idea on counterflow superconductivity with refer- 
ence to electron-hole bilayers was put forward in [10,11]. 
The attempts to observe counterflow conductivity directly 
were done [12-14] for bilayer quantum Hall systems 
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realized in GaAs heterostructures. In the latter systems 
superconducting behavior might be accounted for magne- 
toexcitons [15,16]. The effect is expected for the filling fac- 
tors of Landau levels v; = 7jil 1 ni{l = y/hc/eB is magnetic 
length, «, is the electron density in the ith layer) satisfying 
the condition Vi + v 2 = 1. The role of holes is played by 
empty states in zero Landau level. In experiments [12-14], 
an exponential increase of the counterflow conductivity 
under lowering of temperature was observed, but zero- 
resistance state was not achieved. The latter can be 
explained by the presence of unbound vortices [17-19]. 
Such vortices may appear due to spatial variation of the 
electron density caused by disorder. 

To demonstrate counterflow superconductivity quan- 
tum Hall bilayers should have the parameters that 
satisfy two additional conditions: d < I and I < a* B , 
where d is the interlayer distance, and a* B = e/? 2 /e 2 m* is 
the effective Bohr radius (e is the dielectric constant of 
the matrix, and m* is the effective electron mass). The 
first inequality comes from the dynamical stability con- 
dition. For balanced bilayers (v! = v 2 ) the mean-fields 
theory yields d < 1.175 t. The second inequality is the 
condition for the Coulomb energy e 2 /e£ be smaller than 
the energy distance between Landau levels. In GaAs 
a* B f» 10 nm and the condition I < a B is fulfilled at 
rather strong magnetic fields B> 6T (actually, the 
experiments [12-14] were done at smaller fields). At 
d < 10 nm the interlayer tunneling is not negligible 
small and may result in a locking of the bilayer for the 
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counterfiow transport at small input current [20,21]. At 
larger input current the system unlocks, but the state 
becomes nonstationary one [22-24] that is accompanied 
by a dissipation (the power of losses is proportional to 
the square of the amplitude of the interlayer tunneling 
[22,24]). 

The idea to use graphene for the realization of elec- 
tron-hole superfluidity in quantum Hall bilayers [6-9] 
looks very attractive. The distance between Landau 
Levels in monolayer graphene is proportional to the 
inverse magnetic length, magnetic field does not enter 
into the condition of smallness of the Coulomb energy, 
and small magnetic fields can be used. Smaller magnetic 
fields correspond to smaller critical temperature, but, at 
the same time, they correspond to larger critical d. Use 
of large d allows to suppress completely negative effects 
caused by interlayer tunneling. 

In this article, we concentrated on three questions. 
First, we determine, in what range of internal para- 
meters and external fields magnetoexciton superfluidity 
can be realized. Second, we evaluate critical temperature 
for pure system. Third, we consider its reduction caused 
by electron-impurity interaction. Our study extends the 
results of [8], where a system of two graphene layers 
embedded into a bulk dielectric matrix was considered. 
Here we investigate structures with one and two gra- 
phene layers situated at the surface. 

2 Conditions for the electron-hole pairing in zero 
Landau level 

Quantum Hall effect in graphene is characterized by 
unusual systematics of Landau levels and the additional 
four-fold degeneracy connected with two valleys and 
two spin projections [25]. The energies of Landau levels 

in graphene are £ ±N = ±^^2 |N|, where N = 0, 1, 2, 

and v F ~ 10 6 m/s is the Fermi velocity. In a free 
standing graphene, the N = 0 Landau level is half-filled. 
A state with only completely filled Landau levels corre- 
sponds to a plateau at the Hall conductivity plot (depen- 
dence of a xy on electron density). A free standing 
graphene is just between two plateaus [26]. A given 
quantum states in zero Landau level is characterized by 
the guiding center index X and the combination of the 
spin and valley indexes. Below we call four possible 
combinations, the components, and numerate them by 
the index ji = 1, 2, 3, 4. 

We describe electron-hole pairing in zero Landau level 
in graphene by the wave function that is a generalization 
of the wave function [15] to the multicomponent case 



Here <™ x is the electron creation operator (the opera- 
tor that fills a given state in N = 0 Landau Level), |0) is 
the state with empty zero level, i is the layer index. The 
u - v coefficients satisfy the condition + \vp\ 2 = 0. 
The function (1) can be rewritten in the form 



i*) =nn(^ 

P x 



(2) 



where h\p X 



Ap=UpVp=^l 
mally imbalanced {vp 



C\px is the hole creation operator, and 
the vacuum state is defined as |wfcl) = \\p Tlx c ipx 1 0). 
One can see that the function (2) is an analog of the 
BCS function in the Bardin-Cooper-Schrieffer theory of 
superconductivity. 

The quantity vp = \up \ 2 - \vp \ 2 gives the filling 
factor imbalance for the component [}. The order para- 
meter of the electron-hole pairing reads as 

i>pe" p /2- If a given component is maxi- 

±1) the order parameter Ap is 
equal to zero. 

If a one component bilayer system is balanced, the 
order parameter for the electron-hole pairing is maxi- 
mum. But if the number of components is even, the bal- 
ance v p = 0 can be reached at vp = 1 for half of the 
components and vp = — 1 for the other half. In the latter 
case all Ap = 0. As is shown below, just such a state cor- 
responds to the energy minimum. In other words, in 
balanced graphene bilayers electron-hole pairing does 
not occur. 

At nonzero imbalance ~^2p vp ±2, ±4 at least for 
one component vp =/±l, and electron-hole pairing may 
occur. Nonzero imbalance can be provided by electrical 
field directed perpendicular to the layers. Such a field 
can be created by a voltage difference applied between 
top and bottom gates (see, Figure 1). 

We consider the general structure "dielectric 1-gra- 
phene 1 -dielectric 2-graphene 2-dielectric 3" with three 
different dielectric constants s lt e 2 , and e 3 . Dielectrics 1 
and 3 are assumed to be thick (much thicker than the 
distance between graphene layers d). Solving the stan- 
dard electrostatic problem we obtain the Fourier com- 
ponents of the Coulomb interaction Vw for the 
electrons located in i and i' graphene layers 
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Figure 1 Schematic view of the system under study. C1-C4 are the contacts 



For electrons in N = 0 Landau level in graphene the 
Hamiltonian of Coulomb interaction has the form 



„2»2 



(6) 



X,X' p,p' q 

C lrJX + q r l 2 /2 C l'p'X'~q r i 2 /2 C i'P'^+<1 r l 2 /2 c ipX-q r ey2, 

where S is the area of the system. The interaction with 
the gate field is described by the Hamiltonian 



H G = -— - ^2 (c + lffx c lff x - C2px c ipx), 
xp 



(7) 



where V g is the interlayer voltage created by the exter- 
nal gate (bare voltage). 
Rewriting the wave function (1) in the form 

i*> = n n ( c ° s + ^ sin y c w) i °> - (8) 

x p 

and computing the energy in the state (8) we obtain 

- Jo «>s 2 e/i - ( 2eV s + U T, <™ 0 A , (9) 



cos y« cos b 



where W = e 2 dlF.2l 2 is the energy of direct Coulomb 
interaction. The exchange interaction energies 



determine the parameters / 0 = (Ju + /22)/2 - J12 an d Jz 
= In - Iin- The relation between Op and vp is given by 
equation vp = cos Bp. 



Taking into account the inequalities W >/ 0 , and J n , J22 >In 
(that can be checked direcdy) we find that at V g = 0 the mini- 
mum of (9) is reached at vi = i>2 = 1, V3 = v 4 = —1. It indi- 
cates the absence of electron-hole pairing in balanced 
systems. 

If V g * 0 and belongs to one of the intervals 



nW + J22 — J12 < eV g < (n + 2)W — Jn + Ju, 



(10) 



where n = -4, -2, 0, 2, the energy minimum is reached 
at Vp a ^±1 for one of the components. We will call 
such a component the active one. 

Let us, for instance, consider the interval (10) with n = 
0. Then the energy minimum is reached at 



eV g +J z /2-W 
W-J 0 



The case Vp a = 0 (with maximum order parameter) 
corresponds to the voltage 



eV, 



2 



(11) 



Equation (11) determines the relation between magnetic 
field and the gate voltage V g . To keep vp a = 0 the gate vol- 
tage should be varied synchronically with B. In particular, at 
/z = 0 0?i = £3) the quantities V g and B are linearly related: 



V„ = 



adBc 

S2 



(12) 



where a « 1/137 is the fine structure constant (the 
relation (12) is given in SI units). 

If only the gate voltage or magnetic field is varied, the 
order parameter (and the critical temperature) changes 
nonmonotonically reaching the maximum at the point 
determined by (11). 
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3 Collective mode spectrum and phase diagram 

The components that belong completely to one layer do 
not take part in the pairing. In what follows we consider 
the dynamics of only the active component. 
We describe the active component by the wave function 



*> =n( cos y c ix + 



. » =/2 + e ' iQ ' X+ * X) Sin | C 2,X-Q rP /l) 1 0) (13) 



(here and below we omit the component index). Equa- 
tion (13) describes the state with nonzero counterflow 
currents. To illustrate this statement we neglect for a 
moment the order parameter fluctuations 
{y x = 0,6 x = 6 a ). 

The order parameter is determined by the equation 



A(r) = ^^Wfe(r)(*|cl, x c 2 ,x'|*). 

X,X' 



(14) 



where 



toft 



1 



x r {x - xy 



■e i 2 e 2l 2 



is the single-particle wave function in the coordinate 
representation, L y is the width of the system. 
Substitution (13) into (14) yields 



smy„ -— 
A(r) = e 2 



(15) 



One can see from equation (15) that Q = (Q x , Q y ) is 
the gradient of the phase of the order parameter. 

Computing the energy in the state (13) and neglecting 
the fluctuations we obtain 



8jt£- 



([W - F s {0)]cos 2 0 a - F D {Q)sm 2 0 a ) , (16) 



where 



OO UK. 

FsM = ^ J ph{pqt 2 )\Vn{p) + V 22 (p)]i ~ d&7) 



and 



F D (cj) = J pJ 0 (pql 2 )V 12 {p)e 2 dp. (18) 

o 

Electrical currents can be found from a variation of 
the energy caused by a variation of the vector-potential 



SE 



j,-5Af. 



(19) 



Here A, is the in-plane component of the vector-potential 
in the layer i. To obtain the explicit expression for the varia- 
tion (19) we replace the phase gradient in (16) with the 

e 

gauge-invariant expression Q — — (A p / i — Ap^), where 

Apjr; is the parallel to the graphene layers component of the 
vector potential in the layer i. Then, using (19) one finds the 
currents 



h h h 8nl 2 dQ ' 
At small gradients Q£ « 1 equation (20) is reduced to 

ji = ^PsoQ, (21) 

where coefficient of proportionality between the cur- 
rent and the phase gradient 



esm 2 6 a dF D {Q) 



(20) 



PsO 



32jt- 



-sin 



oo p 2 £ 2 

%j p 3 V 12 (p)e~~dp 



(22) 



is called the zero temperature superfluid stiffness (the 
definition is given in the following section). Since we neglect 
fluctuations, the expression (20) yields the current at T = 0. 

Implying the fluctuations of the amplitude and the 
phase of the order parameter are small one can present 
the energy as 

E = E 0 +E 2 + ... (23) 
The quadratic in fluctuations term can be diagonalized: 

Ei - £ [»>«(-<!)MfK(<r) + \<p{-i)k,M<p{i) - i('»"z(-"/)M<fM</) + «•)]< (24) 
where 



2V S 



<p{4) 



(25) 



2jt£ 2 



-iqX 



are the Fourier components of the fluctuations. 

Equation (24) yields the energy of fluctuations with the 
wave vector directed along the x axis. The component of 
the matrix K can be presented in form independent of 
the choice of the direction of the coordinate axes 

J^(q,Q) = H(q,Q) - F s (|q|) + F D (|Q|) + E{q,Q)cot 2 0 a , (26) 



K„{q,Q) = sin 2 0 £1 S(q,Q), 



(27) 



M<i<Q) = - cos ^ [ p o (lq + Ql) - Fd (lq - QD] /2, (28) 
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where 



In equation (34) 



H(q,Q) = 



i r Mq) + 

2ni 2 L 2 



q 2 l 2 



V 12 (q)cos(|qx Q| I 2 ) \e 



Fd(IQI)- 



F D (|q + Q|) + F D (|q-Q|) 



(29) 



(30) 



The quantities K a p{q) in (24) are expressed in terms of 
(26) as K ap (q) = K a p[q, Q) \ q=qix .. 

The quantity h cos 8 x/2 can be treated as a z-compo- 
nent of the pseudospin and it is canonically conjugated 
with the phase cf> x . The Fourier transformed quantities 
(25) are defined as canonical variables as well. The equa- 
tions of motion for the quantities m z {q) and <fi{q) read as 



dt 

dm z {q) 
dt 



IK^m^q) - iK z<p (q)<p(q), (31) 



"2 K w(<7M<?) - iK^m^q). (32) 
Equation (31) yield the collective mode spectrum 
KWO01QzM +K Z(p (q). Rotating the axes 



(33) 



one obtains the excitation spectrum at general q 



ff (q, Q) = JK„(q, Q)K a {q, Q)+K^(q, Q) 



At Q = 0 the spectrum (33) is isotropic. It can be pre- 
sented in the Bogolyubov form 



^oM = Js q {s q + ye,). 



(34) 



F D {0)-F D {q) 



(35) 



is the kinetic energy (e q ~hq I2M at ql « 1, where M 
is the magnetoexciton mass, see, for instance [27]), and 



Yq = [H(q, 0) - F s {q) + F D {q)]sm 2 0 o 



(36) 



has the sense of the exciton-exciton interaction energy 
(that includes the direct and exchange parts). 

The condition for the dynamical stability of the state 
(13) is the real valueness of the excitation spectrum 
(34). This condition determines the diapason of dli and 
£; where superfluid magnetoexciton state can be realized. 
To be more concrete we consider three types of hetero- 
structures. Type A is a graphene-dielectric-graphene 
sandwich with two graphene layers at the surface, Type 
B is a graphene-dielectric-graphene-dielectric structure 
with one such a layer, and Type C is a system of two 
graphene layers embedded in a dielectric matrix 
(Figure 2). For simplicity, we imply the same dielectric 
constants e for the interfacial layer and the substrate. 

The dynamical stability condition is fulfilled at 
0 < d/i < d c (s), where d c {e) depends on the imbalance 
parameter Vp a = v a . The dependence d c (s) at v a = 0 is 
shown in Figure 3. 

The requirement for the Coulomb energy be smaller 
than the distance between Landau levels yields the restric- 
tion on e. Since we study the pairing in N = 0 Landau level 
we compare the Coulomb energy with the energy distance 
between N = 0 and N = 1 levels co c = s/lkvp/l 
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Figure 3 Phase diagram at v a = 0 for the graphene bilayers of 
A, B, and C type. Solid curves, d c (e); dashed curves, £ c [d/l). 



We have four parameters that characterize the Cou- 
lomb energy W, hi, hi, and hi- At d/£ < d c the largest 
of them is hi (the intralayer exchange interaction in the 
graphene layer at the surface). Therefore, it is natural to 
consider the condition 



hi < u c 



(37) 



as the additional restriction on the parameters. Equa- 
tion (37) can be rewritten as e >s c (d/i). The quantity e c 
can be understood as a critical dielectric constant. The 
dependence e c (d/t) is also shown in Figure 3. 

Two conditions d/l < d c [s) an d £ >e c (dll) determine 
the range of parameters where one can expect a realiza- 
tion of electron-hole pairing and magnetoexciton super- 
fluiduty in graphene bilayer systems. 

4 Critical temperature 

In a bilayer graphene heterostructure with a fixed d the 
magnetoexciton superfluidity can be realized in a wide 
range of magnetic field. Variation of B at fixed gate vol- 
tage results in a change of imbalance of the active com- 
ponent. Simultaneous tuning of V g allows to keep zero 
imbalance v a = 0 and maximum order parameter under 
variation of B. In this section, we study the dependence 
of critical temperature on magnetic field implying such 
a simultaneous tuning. 

Superfluid transition temperature is given by the Bere- 
zinskii-Kostelitz-Thouless equation [15] 



T c = -Ps{T c ), 



(38) 



where p s {T) is the superfluid stiffness at finite tem- 
perature. The superfluid stiffness is defined as the coeffi- 
cient in the expansion of the free energy in the phase 
gradient F = F 0 + f d 2 rp s (yip) 2 /2. In a weakly nonideal 
Bose gas it is equal to p s = h 2 n s /m, where n s is the 
superfluid density. As was shown in previous section, 
superfluid stiffness determines also the supercurrent. 

Taking into account linear excitations we present the 
free energy F = E 0 - TS in the following form 



F = £ 0 + T^ln 



V 



-£2{q, Q)\ 



(39) 



/ 



Expansion of equation (39) yields the following 
expression for the superfluid stiffness 



l v ^r/i( ! H(q 1 Q)\| 



;N,(1+N,) 



]■ 

' Iq=oJ 



(40) 



It follows from (40) and (33) that p s (T) <p s0 (thermal 
fluctuations reduce the superfluid stiffness). 

For the spectrum I2(q) = E(q) + hqv (where v = tN(j)l 
m is the superfluid velocity) (40) yields the well-known 
answer for the superfluid density [28]. Equation (40) 
generalizes the results [28] for the general case. 

The dependence of critical temperature on magnetic 
field at v a = 0 and e = 4 is shown in Figure 4. One can 
see that the maximum critical temperature is reached 
approximately at B « 0.5-B^, where Bj = q>/nd 2 with <p = 
hc/2e, the magnetic flux quantum. 

5 Influence of impurities on the critical 
parameters 

In the previous section, we have determined the influ- 
ence of thermal fluctuations on the superfluid stiffness. 
In this section, we consider the effect of reduction of 




B/B d 

Figure 4 Critical temperature vs magnetic field for A, B, and C 
structures. Temperature is given in units of T d = e 2 /£<J, 
magnetic field, in units of B d = cp/nd 2 
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the superfluid stiffness caused by the interaction of mag- 
netoexcitons with impurities. 

The Hamiltonian of the interaction of the active com- 
ponent with impurities can be presented in the form 



(41) 



where U z (q) = £/i(q) - U 2 (q), Uj(q) is the Fourier- 
component of the impurity potential in the layer i, and 

/ ^2^2 \ 

Mq) = E c+ qy e c ci Y e^{-^ x - — k) 

X iv.lL iY—— v ' 



2 2 



is the Fourier component of the electron density 
operator for the active component. 

In the state (13), the energy of interaction with the 
impurities expressed in terms of m z {q) reads as 



const + ^U z [q)m z {q), 

1 

where 



(43) 



^2jt£ 2 S 



The interaction (43) induces the fluctuations of the 
density and the phase of the order parameter. 

Their values can be obtained from the Euler-Lagrange 
equations 



SE 



Sm z (q) 
SE 

M4) 



= 0, 



0, 



(44) 



where E is the energy of the system, described by the 
Hamiltonian H = H c + H G + H^p in the state (13). 

Equations (44) solved in linear in impurity potential 
approximation yield 



m z {q) 



1 K„{q)U z {qk) 



iK Zf {q)U z {(pL) 



(45) 



(46) 



Substituting (45), (46) into the expression for the 
energy one finds the correction to the energy caused by 
the electron-impurity interaction 



In equation (47), the contribution of fluctuations with 
the wave vectors directed along x is taken into account. 
Summing the contribution for all wave vectors one obtains 



qH 1 

Ab = 83r£ 2 S Z^™ q V( q ,Q)K w ( q ,Q)-/^(q,Q)' 



(48) 



For simplicity, we specify the case where impurities 
are located in graphene layers. Then the Fourier-compo- 
nent of the impurity potential can be presented in the 
form 



(49) 



where t a are the impurity coordinates, and u z j(q) = U\ t 
,(q) - u 2 ,i{q) with M^,-(q), the potential in the layer k of a 
single impurity centered at r = 0 in the layer i. 

Averaging over impurities yields 

qV 

K„{q,Q)e~ 2 (50) 



Q) 



where « imp is the impurity concentration in a layer. 
At Q£ « 1 the energy (50) can be expanded in series 



AE = AE Q + -Ap s Q 2 , 



where 



k.Mr + Mi)! 



K'Jq.O) 



(51) 



(52) 



is the correction of the superfluid stiffness. One can 
check that the correction Ap s is negative. Thus, the 
interaction with impurities results in decrease of critical 
parameters. 

At v a = 0 equation (52) is reduced to 



qV 

,(ki(q)| 2 -<-Mq)| 2 )^~ 



qV 



(53) 



where p s0 (equation (22)) is taken at 8 a = nil. 

The shift of critical temperature is evaluated as A TJ T c 
~ Ap s /p s0 .* We define the critical impurity concentration 
n i m p as a concentration at which Ap s /p s0 = 1. We con- 
sider charged impurities with the potential u^iq) = (-1)' 
(^i2(q) - V,'i(q))- The dependence of critical impurity 
concentration on magnetic field at e = 4 and v a = 0 is 
shown in Figure 5. We also evaluated critical concentra- 
tions for neutral impurities. These concentrations are 
much larger, and the influence of neutral impurities can 
be neglected. 
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B/B d 

Figure 5 Critical impurity concentration versus magnetic field 
for charged impurities located in graphene layers 



6 Conclusion 

In conclusion, we present some estimates. Let us specify 
the type B structure (the one used in [1]) with d = 20 
nm and e = 4. For this structure the maximum critical 
temperature T c ~ 3 K (in pure case) is reached in mag- 
netic field B ~ 0.8 T. At such B the critical impurity 
concentration is n\ mp « 2 • 10 9 cm~ 2 . The gate voltage 
determined by equation (11) is V g ~ 6 mV, that corre- 
sponds to electrostatic field E ~ 3 kVcm 1 . 

Basing on the results of our study we may state the 
following. 

1. Graphene bilayer structures are perspective 
objects for the observation of magnetoexciton super- 
fluidity. The advantages are smaller magnetic fields 
and no restriction from above on physical interlayer 
distance, that means the possibility to suppress com- 
pletely interlayer tunneling. 

2. Gate voltage should be created between graphene 
layers for a realization of magnetoexciton 
superfluidity. 

3. Certain conditions on dielectric constant and on 
the ratio between interlayer distance and magnetic 
length should be satisfied. 

4. Structures with graphene layers situated at the 
surface have larger critical parameters. 

5. Neutral impurities are not dangerous for the mag- 
netoexciton superfluidity, but the concentration of 
charged impurities should be controlled. 

Endnote 

a Since in our approach we assume smallness of Ap s lp s0 
it is just an estimate. 
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